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A  kinetic  model  to  predict  the  relative  intensities  of  the  atomic  C/H/N/O  emission  lines  in  laser-induced 
breakdown  spectroscopy  (LIBS)  has  been  developed  for  organic  compounds.  The  model  includes  a  com¬ 
prehensive  set  of  chemical  processes  involving  both  neutral  and  ionic  chemistry  and  physical  excitation 
and  de-excitation  of  atomic  levels  affecting  the  neutral,  ionic,  and  excited-state  species  concentrations. 
The  relative  excited-state  atom  concentrations  predicted  by  this  modeling  are  compared  with  those  de¬ 
rived  from  the  observed  LIBS  intensities  for  355  nm  ns  laser  irradiation  of  residues  of  two  organic  com¬ 
pounds  on  aluminum  substrate.  The  model  reasonably  predicts  the  relative  excited-state  concentrations, 
as  well  as  their  time  profiles.  Comparison  of  measured  and  computed  concentrations  has  also  allowed  an 
estimation  of  the  degree  of  air  entrainment.  ©  2010  Optical  Society  of  America 
OCIS  codes:  300.6365,  000.1570,  350.5400. 


1.  Introduction 

There  has  been  considerable  interest  in  the  emission 
spectra  from  laser-induced  plasmas.  From  an  analyti¬ 
cal  point  of  view,  the  relative  intensities  of  the  various 
atomic  and  molecular  features  provide  information 
with  which  to  identify  the  material  from  which  the 
plasma  was  created.  The  techniques  of  laser-induced 
breakdown  spectroscopy  (LIBS)  [1-4]  and  laser  abla¬ 
tion  coupled  plasma  mass  spectrometry  [5]  are  based 
on  this  idea.  Because  little  or  no  sample  preparation  is 
required,  LIBS  has  been  developed  into  a  powerful 
method  for  the  identification  of  many  materials. 

The  relative  intensities  of  the  various  atomic  and 
molecular  spectroscopic  signatures  depend  upon 
the  material  entrained  into  the  plasma  and  the  phy¬ 
sical  and  chemical  processes  occurring  within  the 
plasma.  We  are  interested  in  developing  a  full  chemi¬ 
cal  and  physical  model  of  the  LIBS  plume  for  various 
systems  in  order  to  understand  in  detail  the  various 
processes  leading  to  the  observable  spectroscopic  sig- 
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natures.  A  full  treatment  of  the  processes  occurring 
within  the  plume  involves  an  understanding  of  both 
the  plasma  hydrodynamics  and  the  relevant  chemical 
and  physical  processes.  In  this  work  on  LIBS  of  two 
representative  organic  compounds,  we  concentrate 
on  the  latter  processes  and  use  experimental  observa¬ 
tions  of  the  plasma  to  constrain  its  temperature  and 
electron  density.  We  present  a  kinetic  model  to  de¬ 
scribe  the  neutral  and  ionic  chemistry,  as  well  as  pro¬ 
cesses  leading  to  the  production  and  removal  of  the 
excited,  emitting  C/H/N/O  atomic  levels.  In  previous 
work,  we  carried  out  similar  studies  of  LIBS  of 
metallic  lead  in  air  and  argon  atmospheres  [6,7],  as 
well  as  LIBS  of  RDX  in  air  [8].  In  the  present  study, 
we  compare  the  results  of  our  kinetic  modeling  with 
experimentally  determined  relative  LIBS  emission 
intensities  in  order  to  test  the  ability  of  our  model 
to  describe  the  excited-state  populations. 

Several  groups  have  developed  models  to  explain 
the  plasma  formation  and  plume  expansion  dynamics, 
and  representative  studies  are  described  here.  Bogae- 
rts  and  co-workers  have  considered  the  laser  ablation 
of  Cu  and  the  plasma  expansion  in  different  back¬ 
ground  gases  and  in  the  double-pulse  configuration 
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[9,10].  Capitelli  and  co-workers  have  considered 
both  equilibrium  and  nonequilibrium  conditions  in 
laser-induced  plasmas  through  a  time-dependent 
collisional-radiative  model  of  the  atomic  excited-state 
populations  and  the  electron  energy  distribution 
[11,12] .  Vidal  et  al.  [13]  carried  out  a  one-dimensional 
Lagrangian  hydrodynamics  calculation  of  the  plasma 
expansion  from  ablation  of  aluminum  in  ambient 
air.  Ashfold  and  co-workers  have  carried  out  2- 
dimensional  plasma  and  chemical  modeling  of  a  dia¬ 
mond  deposition  reactor  and  compared  the  results 
with  experimental  observations  [14]. 

2.  Kinetic  Model 

We  have  considerably  extended  the  kinetic  model  pre¬ 
viously  developed  [6]  to  describe  air  oxidation  and  io¬ 
nization  in  the  LIBS  plume  of  lead  in  air.  Our  kinetic 
model  includes  reactions  describing  air  oxidation  and 
plasma  chemistry,  as  well  as  the  thermal  decomposi¬ 
tion  of  small  organic  compounds  (those  containing 
<4  carbon  atoms).  The  rate  constants  were  obtained 
from  various  databases  and  review  articles  [15-20], 
as  well  as  a  number  of  research  papers. 

A  significant  aspect  of  our  kinetic  model  is  the  in¬ 
clusion  of  the  excited  electronic  states  of  H,  C,  N,  and 
0  atoms,  which  are  required  to  describe  the  popula¬ 
tions  of  the  emitting,  excited  states  of  these  ele¬ 
ments.  Table  1  lists  the  atomic  states  included  in  the 
model,  as  well  as  the  relevant  emission  lines.  We 
have  included  transitions  in  the  VUV,  which  are 
not  observed  in  LIBS  spectra  but  which  depopulate 
the  lower  levels  of  the  observed  LIBS  lines.  While  not 
radiatively  connected  with  the  observed  emission 
lines,  the  low-lying  C,  N,  and  0  levels  associated  with 
the  ground  electron  configurations  are  included  since 
their  concentrations  can  be  significant  and  electron- 
impact  excitation  of  these  levels  can  populate  the 
emitting  levels. 

We  now  consider  the  process  populating  and  de¬ 
pleting  the  excited  atomic  electronic  states.  The  ex¬ 
cited  states  were  assumed  to  be  populated  only  by 
electron-impact  excitation  of  the  states  associated 
with  the  ground  atomic  electron  configuration.  The 
temperature-dependent  excitation  rate  constants  for 
H  and  C  were  taken  from  the  National  Institute  for 
Fusion  Science  atomic  and  molecular  database  [21], 
while  the  corresponding  rate  constants  for  N  and  0 
were  derived  from  effective  collision  strengths  com¬ 
puted  by  Tayal  [22]  and  Zatsarinny  and  Tayal  [23], 
respectively.  Electron-impact  de-excitation  rate  con¬ 
stants  were  computed  by  detailed  balance.  Heavy- 
particle  excitation  was  not  included  in  the  model 
since  this  process  occurs  predominantly  at  keV  ener¬ 
gies  [24];  similarly,  production  of  excited  states  by 
electron-ion  recombination  was  neglected.  In  addi¬ 
tion  to  electron-impact  de-excitation,  radiative  decay 
[25]  and  collisional  electronic  quenching  by  neutral 
species  were  included  as  processes  depopulating 
excited  atomic  states.  Unfortunately,  electronic  qu¬ 
enching  rate  constants  are  available  only  for  stable 
collision  partners  at  room  temperature. 


Table  1.  Atomic  Energy  Levels  Included  in  the  Kinetic  Model 


Notation 

Level 

Energy  (cm  4) 

Emission  Line  (nm)“ 

Hydrogen 

H 

n  —  1 

0 

97.2,  102.6,  121.6  (L) 

H*2 

n  =  2 

82259 

121.6  (U),  486,  656  (L) 

H*3 

n  =  3 

97492 

102.6  (U),  656  (U) 

H*4 

n  =  4 

102824 

97.2,  486  (U) 

Carbon 

C 

2s2  2p2  3P 

0 

156.1,  165.7  (L) 

C*2 

2s22pz  1D 

10194 

193.7  (L) 

C*3 

2s2  2p2  4S 

21648 

247.9  (L) 

C*4 

2s22p3s  3P° 

60360 

165.7  (U) 

C*5 

2s22p3s  1P° 

61982 

193.7,  247.9  (U) 

C*6 

2s2p3  3D° 

64091 

156.1  (U) 

Nitrogen 

N 

2 p3  4S° 

0 

120.0,  113.4  (L) 

N(2D) 

2 p3  2D° 

19226 

149.3  (L) 

N(2P) 

2 p3  2P° 

28840 

174.3  (L) 

N*4 

2p23s  4P 

83337 

120.0  (U),  (742-746), 

N*5 

2pz3s  2P 

86192 

(818-821), (868-874)  (L) 
149.3,174.3  (U) 

N*6 

2s2p4 4P 

88134 

113.4  (U) 

N*86 

2p23p  4D° 

94839 

(868-874) (U) 

N*9 

2p23p  4P° 

95511 

(818-821)  (U) 

N*10 

2p23p  4S° 

96752 

(742-746) (U) 

Oxygen 

0 

2 p4  3P 

0 

130.4  (L) 

0(1D) 

2 p4  1D 

15868 

0(1S) 

2 p4  LS 

33793 

0*4 

3s5S° 

73768 

777  (L) 

0*5 

3s3S° 

76794 

130.4  (U),  845  (L) 

0*6 

3 p5P 

86626 

777  (U) 

0*7 

3 p3P 

88631 

845  (U) 

“U  or  L  denotes  that  this  level  is  the  upper  or  lower  level,  respec¬ 
tively,  of  the  emission  line. 

6N*7  denotes  the  2p23p2S°  state,  which  lies  at  93582  cm-1  but 
which  was  not  included  in  the  kinetic  model. 

While  CN  and  C2  molecular  emission  is  usually 
observed  in  LIBS  of  organic  compounds  [26,27],  as 
is  the  case  in  the  present  experiments,  the  available 
data  on  electron-impact  excitation  of  the  excited  mo¬ 
lecular  states  are  lacking,  and  we  did  not  include 
these  states  in  our  model.  Since  our  organic  com¬ 
pounds  were  prepared  as  residues  on  Al  foil  (see  Sec¬ 
tion  3),  emission  of  Al  lines,  due  to  ablation  of  the 
substrate,  was  also  observed.  Sufficient  literature 
data  are  available  to  include  several  Al  species  in 
our  model  [Al,  ALAS,  Al+,  AlO,  A1H]  in  order  to  allow 
us  to  interpret  the  Al  394-396  nm  [4s2S]/2  -> 
3p2P°1/2,3/2]  line  intensities  relative  to  those  of  the 
light  atoms.  The  energy-dependent  electron-impact 
excitation  cross  section  to  the  4S  state  was  taken  from 
Smirnov  [28].  The  early  stage  of  Al  oxidation  was 
taken  into  account  by  inclusion  of  the  species  AlO 
and  A1H. 

The  Chemkin  package  [29]  was  employed  to  solve 
the  differential  equations  governing  the  time- 
dependent  concentrations.  We  have  ignored  spatial 
variations  of  the  concentrations  and  have  assumed 
a  homogeneous  closed  reactor  model.  To  match  the 
required  format  of  the  input  data,  the  temperature 
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dependence  of  the  rate  constants  of  the  included 
processes  were  fitted  to  the  functional  form  k  =  AT" 
exp (~Ea/RT),  where  T  is  the  temperature  (in  de¬ 
grees  Kelvin)  and  A,  n,  and  E„  are  fitted  parameters. 
A  total  of  134  species  (including  the  atomic  states 
listed  in  Table  1)  and  1501  processes  (with  forward 
and  reverse  processes  listed  separately)  were  in¬ 
cluded  in  the  kinetic  model.  A  list  of  the  processes 
and  the  kinetic  parameters,  as  well  as  the  sources 
for  the  rate  constants,  can  be  obtained  from  one  of 
the  authors  (PJD).  In  addition  to  the  matrix  of  rate 
constants,  integration  of  the  rate  constants  requires 
specification  of  the  initial  concentrations,  which  are 
determined  by  the  laser-matter  interaction,  and 
the  temperature  temporal  profile.  As  described  in 
Section  4,  we  make  various  assumptions  about  the 
initial  concentrations  and  derive  an  estimated  tem¬ 
perature  profile  from  experimental  data.  In  all  cases, 
we  assumed  a  constant  pressure  of  latm  in  the 
simulations. 

3.  Experimental  Apparatus 

The  apparatus  employed  for  measurement  of  LIBS 
intensities  is  similar  to  that  used  in  previous  LIBS 
investigations  in  our  laboratory  [27,30],  and  a  brief 
description  is  provided  here.  Single  pulses  of  355  nm 
radiation  (tripled  Nd:YAG  fundamental  from  a 
Continuum  Powerlite  8000,  ~7  ns  pulse  width,  10  Hz 
repetition  rate)  were  passed  through  a  mechanical 
shutter  and  focused  with  a  50  mm  focal  length  lens 
onto  the  sample,  which  was  mounted  on  a  motorized 
translation  stage.  A  fresh  area  was  irradiated  on 
each  laser  shot  in  ambient  air.  Emission  (200- 
860  nm)  from  the  LIBS  plasma  was  focused  with  a 
pair  of  30°  off-axis  parabolic  mirrors  (focal  length 
25.4  mm)  onto  a  UV-transmitting  optical  fiber  (50 /mi 
core).  The  output  of  the  fiber  was  transmitted  to 
an  echelle  spectrometer  (Andor  Mechelle  5000) 
equipped  with  a  gated,  intensified  charge  coupled  de¬ 
vice  (ICCD)  (Andor  DH734- 18-03).  The  ICCD  delay 
and  gate  width  were  controlled  by  a  digital  delay  gen¬ 
erator  (Stanford  Research  Systems  DG535).  The  gate 
width  was  fixed  at  100  ns,  and  delays  from  20  ns  to 
3  /is  were  employed.  The  wavelength-dependent  sen¬ 
sitivity  of  the  optical  fiber/spectrometer/ICCD  combi¬ 
nation  was  determined  by  recording  spectra  of 
calibrated  deuterium  and  tungsten-halogen  lamps, 
and  the  spectra  were  corrected  for  the  varying  sensi¬ 
tivity.  Atomic  line  intensities  were  computed  by  de¬ 
termining  the  area  under  the  line  profile  above  any 
background. 

The  two  organic  compounds  investigated,  anthra¬ 
cene  and  2, 4-dinitro toluene  (Sigma-Aldrich),  were 
prepared  as  residues  on  aluminum  foil  substrates. 
For  each  compound,  a  10 /d  aliquot  of  a  nearly 
saturated  solution  in  methanol  was  delivered  to 
the  substrate  surface,  and  the  solution  was  allowed 
to  dry  in  the  hood.  The  resulting  spot  size  was  mea¬ 
sured  to  determine  the  average  surface  coverage, 
which  was  ~25  and  40  //g/cm2  for  anthracene  and 
2,4-dinitrotoluene,  respectively. 


4.  Results  and  Discussion 

A.  Bare  Al  Foil 

We  first  consider  LIBS  of  bare  Al  foil  in  ambient  air 
in  order  to  investigate  the  relative  concentrations  of 
the  light  atoms  and  the  Al  atom  in  excited,  emitting 
states.  We  recorded  LIBS  spectra  with  the  detector 
gate  delay  stepped  in  increments  of  100  ns.  The  laser 
pulse  energy  was  20mJ  for  these  measurements. 
Figure  1  displays  the  Al  394-396  nm  lines  at  several 
detector  gate  delays.  It  can  be  seen  that  the  lines  are 
quite  broad  at  short  delays  and  the  widths  decrease 
with  increasing  gate  delay.  The  width  of  the  lines  is  a 
result  of  Stark  broadening  due  to  the  high  initial 
electron  density.  It  can  be  seen  in  Figs.  1(b) — 1(d)  that 
at  longer  delays  there  is  self-absorption  of  the  lines, 
due  to  the  presence  of  Al  atoms  in  the  cooler  parts  of 
the  plasma.  It  is  also  evident  that  the  lines  are  Stark 
shifted,  as  well  as  broadened. 

The  wings  of  the  line  profiles  were  fitted  to  a  Voigt 
function  in  order  to  determine  the  Stark  width  as  a 
function  of  time.  The  fitted  widths  were  employed  to 
estimate  the  electron  density,  using  literature  data 
[31,32]  for  the  linear  Stark  broadening  parameter 
of  the  394  nm  line.  The  derived  electron  density  as 
a  function  of  time  is  presented  in  Fig.  2(b).  It  can 
be  seen  that  the  electron  density  at  the  shortest 
delay  is  estimated  to  be  ~3  x  1018  electrons  cm-3. 
It  should  be  noted  that  this  electron  density  is  ap¬ 
proximately  an  order  of  magnitude  greater  than  the 
density  range  at  which  the  Stark  parameters  were 
measured,  and  hence  our  derived  densities  could 
have  a  systematic  error. 

It  is  also  of  interest  to  estimate  the  plasma  tem¬ 
perature  as  a  function  of  time,  since  this  is  important 
input  for  our  kinetic  modeling,  as  discussed  above. 
We  have  utilized  the  ratio  of  the  intensities  of  the 
Al  394.40  and  308.22  nm  lines  and  the  assumption  of 
local  thermodynamic  equilibrium  to  derive  the  tem¬ 
perature  as  a  function  of  time.  In  this  analysis,  the 
line  intensities  were  computed  as  areas  of  Voigt  pro¬ 
files  fitted  to  the  wings  of  the  lines,  as  both  lines  were 
subject  to  self-absorption  (see  Fig.  1).  Figure  2(a)  pre¬ 
sents  the  derived  temperature  as  a  function  of  time. 
It  can  be  seen  that  the  temperature  at  the  shortest 
delay  is  ~18,000K  and  drops  to  ~10,000K  by  3/^s. 
These  values  were  least  squares  fitted  to  an  exponen¬ 
tial  decay  plus  a  constant  term  in  order  to  have  a 
smooth  temperature  profile  as  input  for  the  kinetic 
simulations.  The  scatter  of  the  points  about  the  fitted 
line  gives  an  indication  of  the  precision  of  this  tem¬ 
perature  measurement.  It  is  certainly  likely  that  the 
plasma  at  the  shortest  delays  is  not  in  thermody¬ 
namic  equilibrium;  however,  this  is  the  simplest 
model  for  use  for  input  in  the  kinetic  simulations. 

Figure  3  presents  the  observed  time  dependence  of 
the  relative  concentrations  of  the  excited,  emitting 
states  of  the  light  atoms  and  Al.  These  concentra¬ 
tions  were  determined  from  the  areas  of  the  emission 
lines  and  the  known  [25]  radiative  decay  rates;  the 
line  intensity  is  proportional  to  the  excited-state 
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Fig.  1.  Spectra  of  the  A1 394-396  nm  lines  in  LIBS  of  bare  A1  foil  at  detector  gate  delays  of  (a)  20  ns,  (b)  500  ns,  (c)  1000  ns,  and  (d)  2000  ns, 
with  a  gate  width  of  100  ns.  The  dotted  lines  represent  Voigt  profile  fits  to  the  wings  of  the  lines.  The  355  nm  laser  energy  was  20  mj  per 
pulse. 


x  10 


Fig.  2.  (a)  Plasma  temperature  and  (b)  electron  density  as  a  func¬ 
tion  of  time  after  the  355  nm  laser  pulse  (20  mj)  for  LIBS  of  bare  A1 
foil.  The  former  was  derived  from  the  relative  intensities  of  the  A1 
394.40  and  308.22  nm  lines,  while  the  latter  was  obtained  from  the 
Stark  broadening  of  the  394.40  nm  line  (see  text). 


concentration  multiplied  by  the  radiative  decay  rate. 
We  see  that  the  concentrations  of  the  excited  N  and  0 
states  are  comparable  in  magnitude  and  have  a  simi¬ 
lar  time  dependence.  The  initial  concentration  of  the 
A1  AS  state  is  somewhat  greater,  and  the  concentra¬ 
tion  of  this  state  decays  more  slowly  with  time  than 
the  concentrations  of  the  N  and  O  states. 

In  order  to  interpret  the  observed  excited-state  con¬ 
centrations,  we  have  carried  out  kinetic  simulations 
of  the  time-dependent  species  concentrations  with 
methods  described  in  Section  2.  The  temperature 


N*10,  0*6,  and  0*7  levels  (in  the  notation  of  Table  1)  determined 
from  the  measured  intensities  (corrected  for  the  wavelength- 
dependent  detection  sensitivity)  of  atomic  emission  lines  for  the 
355  nm  laser  irradiation  (pulse  energy  20  mj)  of  bare  A1  foil  in  am¬ 
bient  air. 
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profile  in  these  simulations  was  taken  from  the  fit 
[smooth  line  in  Fig.  2(a)]  to  the  temperatures  esti¬ 
mated  from  the  A1  atomic  line  ratios.  As  discussed 
in  Section  2,  one  critical  assumption  in  this  modeling 
is  the  species  concentrations  after  the  laser  pulse.  We 
must  also  make  assumptions  as  to  the  entrainment  of 
the  ablated  A1  and  the  ambient  air  into  the  LIBS 
plume.  We  have  carried  out  simulations  for  various 
A1 :  air  mole  ratios  in  the  plasma  and  find  that  our  ob¬ 
served  ratio  of  the  ALLS’  concentration  to  that  of  the 
excited  N  and  0  states  is  consistent  with  a  very  small 
A1 :  air  mole  ratio  in  the  plasma. 

We  present  in  Fig.  4  computed  time-dependent 
concentrations  of  the  excited,  emitting  states  for 
an  assumed  A1 :  air  mole  ratio  of  0.1  and  various  as¬ 
sumptions  for  the  initial  species  concentrations.  In 
all  cases,  the  computed  AL4S  concentration  is  sev¬ 
eral  orders  of  magnitude  higher  than  those  of  the 
N  and  0  states.  In  large  part,  this  is  due  to  the  much 
larger  rate  constant  for  electron-impact  excitation  of 
the  A1  4S  state  than  those  for  the  excited,  emitting  N 
and  0  states;  the  difference  in  the  magnitudes  of  the 
rate  constants  is  due  to  the  lower  Al  4S  excitation 
energy. 

It  can  be  seen  in  Fig.  4  that  the  excited-state  con¬ 
centrations  reach  their  maximum  values  at  greatly 
different  times  for  the  various  assumed  initial  con¬ 
centrations.  For  assumed  simple  heating  [Fig.  4(a)] 
the  maximum  N  and  0  concentrations  are  reached 
~300ns  after  the  laser  pulse.  In  this  case,  the 
A14S  concentration  is  computed  to  increase  slightly 


after  this  induction  period.  By  contrast,  the  maxi¬ 
mum  concentrations  are  reached  ~  100  ns  after  the 
laser  pulse  when  complete  atomization  is  assumed 
[Fig.  4(b)],  We  see  a  similar  concentration  temporal 
profile  when  10%  ionization  of  A1  is  additionally  as¬ 
sumed;  the  relatively  small  amounts  of  ions  and  elec¬ 
trons  with  this  assumption  do  not  significantly  affect 
the  computed  concentrations.  The  computed  concen¬ 
trations  with  the  assumption  of  complete  atomiza¬ 
tion  and  10%  ionization  of  N  are  seen  in  Fig.  4(d) 
to  lead  to  computed  excited  N  and  O  atom  time- 
dependent  concentrations  in  reasonable  agreement 
with  experiment.  In  this  case,  the  excited-state  con¬ 
centrations  rise  rapidly  after  the  laser  pulse  on  the 
time  scale  of  the  plot  and  decay  with  a  rate  similar 
to  that  seen  in  the  experimentally  determined  ex¬ 
cited-state  concentrations  (Fig.  3).  Additionally,  the 
calculations  reproduce  the  observed  much  slower 
decay  of  the  AT4S  state. 

B.  Anthracene 

We  chose  anthracene  (C14H10)  as  one  of  the  organic 
compounds  for  study  since  the  relative  intensities  of 
the  N  and  O  lines  are  a  good  measure  of  the  degree  of 
air  entrainment  for  this  hydrocarbon.  Figure  5(a) 
presents  the  determined  time  dependence  of  the  re¬ 
lative  concentrations  of  the  excited,  emitting  states 
of  the  light  atoms  and  A1  for  LIBS  of  residues  of  an¬ 
thracene  on  A1  foil.  A  laser  pulse  energy  of  20  mJ  was 
employed.  The  concentrations  and  time  dependence 
of  the  excited  C,  N,  and  O  states  are  seen  to  be  simi¬ 
lar.  The  concentration  of  the  H*3  state  is  similar  to 


Fig.  4.  (Color  online)  Plots  of  the  time-dependent  concentrations  of  the  excited,  emitting  atomic  levels  (in  the  notation  of  Table  1)  for  an 
Al-air  mixture  (A1 :  air  mole  ratio  0.1),  computed  with  the  kinetic  model  described  in  Section  2  and  the  temperature  plotted  in  Fig.  2(a).  The 
assumed  initial  conditions  and  species  mole  ratios  were  as  follows:  (a)  simple  heating  [Al,  1;  N2,  7;  02,  3],  (b)  complete  atomization  [Al,  1;  N, 
14;  O,  6],  (c)  atomization  and  10%  ionization  of  Al  [Al,  0.9;  Al+,  0.1;  E,  0.1;  N,  14;  O,  6],  (d)  atomization  and  10%  ionization  of  N  [Al,  1;N,  12.6; 
N+,  1.4;  E,  1.4;  O,  6].  The  electron  and  A1*4S  concentrations  are  scaled  as  indicated  in  the  panels. 
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Fig.  5.  (Color  online)  (a)  Relative  populations  of  the  excited,  emit¬ 
ting  levels  (in  the  notation  of  Table  1)  determined  from  the  mea¬ 
sured  intensities  (corrected  for  the  wavelength-dependent 
detection  sensitivity)  of  atomic  emission  lines  for  the  355  nm  laser 
irradiation  (pulse  energy  20  mJ)  of  amthracene  residue  on  A1  foil  in 
ambient  air.  (b)  Plasma  temperature,  derived  from  the  relative  in¬ 
tensities  of  the  A1  394.40  and  308.22  nm  lines,  as  a  function  of  time 
after  the  355  nm  laser  pulse  (20  mJ)  for  LIBS  of  anthracene  resi¬ 
due  on  A1  foil. 

that  of  Al*  AS  at  short  delays  but  decreases  more  ra¬ 
pidly.  As  in  the  case  of  LIBS  of  bare  Al  foil,  the  Al  4S 
[Fig.  3],  the  Al  4S  concentration  persists  for  a  signif¬ 
icantly  longer  time  than  do  the  concentrations  of  the 
light  atoms.  We  also  see  that  the  Al  4S  concentration 
extends  to  longer  delays  than  the  excited-state  con¬ 
centrations  of  the  lighter  atoms,  as  was  also  observed 
in  Fig.  3  for  LIBS  of  bare  Al  foil  in  ambient  air. 

Figure  5(b)  displays  the  temporal  dependence  of 
the  plasma  temperature,  estimated  from  the  ratio 
of  the  Al  394.40  and  308.22  nm  lines,  as  determined 
by  the  procedure  described  in  Subsection  4.A.  The  in¬ 
itial  temperature  was  found  to  be  somewhat  higher 
than  in  the  case  of  LIBS  of  bare  Al  foil  [see  Fig.  2(a)]. 

We  have  carried  out  kinetic  simulations  of 
anthracene-air  mixtures  with  the  methods  describ¬ 
ed  in  Section  2.  We  have  not  included  Al  in  the  re¬ 
action  mixture  since  the  Al  mole  fraction  required  to 
model  the  observed  Al'  4S  relative  concentrations 
was  found  in  Subsection  4.A  to  be  much  too  low 
to  affect  the  chemistry  of  the  light  atoms.  The  tem¬ 
perature  profile  in  the  simulations  was  taken  from 
the  fit  [smooth  line  in  Fig.  5(b)]  to  the  temperatures 
estimated  from  the  Al  atomic  line  ratios. 

As  discussed  above,  we  must  also  make  an  assum¬ 
ption  about  the  initial  concentrations,  immediately 
after  the  laser  pulse.  In  the  case  of  LIBS  of  bare 
Al  foil  [Subsection  4. A],  we  found  that  prompt  forma¬ 


tion  of  the  excited,  emitting  states  was  computed 
only  if  complete  atomization  and  partial  ionization 
was  assumed.  We  have  employed  this  assumption 
in  our  kinetic  simulations  for  both  anthracene  and 
2,4-dinitrotoluene  [discussed  in  Subsection  4.C],  In 
the  case  of  anthracene,  we  assume  10%  ionization 
of  the  C  atoms. 

Figure  6  presents  computed  time-dependent  con¬ 
centrations  of  the  excited,  emitting  light  atoms  for 
two  assumed  anthracene-air  mole  ratios  (1:5  and 
1 : 10).  In  both  cases,  the  H*  3  state  is  computed  to  have 
the  highest  concentration  of  the  emitting  states.  The 
second  most  populated  state  is  predicted  to  be  the  C*5 
state.  The  concentrations  of  the  emitting  N  and  O 
states  are  lower  and  are  found  to  scale  approximately 
linearly  with  the  air-anthracene  mole  ratio.  Also 
plotted  in  both  panels  of  Fig.  6  is  the  computed  elec¬ 
tron  density,  which  is  estimated  to  be  several  orders  of 
magnitude  larger  than  the  excited-state  atom  concen¬ 
trations.  We  also  see  that  the  electron  density  is  com¬ 
puted  initially  to  increase  slightly  and  then  decay 
much  more  slowly  than  the  concentrations  of  the 
emitting  states  of  the  light  atoms. 


Fig.  6.  (Color  online)  Plots  of  the  time-dependent  concentrations 
of  the  excited,  emitting  atomic  levels  (in  the  notation  of  Table  1)  for 
an  anthracene— air  mixture,  computed  with  the  kinetic  model  de¬ 
scribed  in  Section  2  and  the  temperature  plotted  in  Fig.  5(a)  .  The 
assumed  initial  conditions  and  species  mole  ratios  were  as  follows: 
(a)  anthracene-air  mixture  at  a  1 : 5  mole  ratio,  with  assumed  com¬ 
plete  atomization  and  10%  ionization  of  C  [C,  12.6;  C+,  1.4;  E,  1.4; 
H,  10;  N,  7;  O,  3]  and  (b)  anthracene-air  mixture  at  a  1 : 10  mole 
ratio,  with  assumed  complete  atomization  and  10%  ionization  of  C 
[C,  12.6;  C+,  1.4;  E,  1.4;  H,10;  N,  14;  O,  6],  The  electron  concentra¬ 
tion  is  scaled  as  indicated  in  the  panels. 
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In  simulations  with  other  temperature  profiles,  we 
found  that  the  relative  concentration  of  the  C*5  state 
is  somewhat  sensitive  to  the  plasma  temperature.  By 
contrast,  the  relative  concentrations  of  the  H,  N,  and 
0  excited  states  do  not  seem  to  depend  significantly 
upon  the  plasma  temperature.  The  sensitivity  of  the 
C*5  state  concentration  is  due  to  the  low  ionization 
potential  of  C  and  the  relatively  low  excitation  en¬ 
ergy  of  the  C*5  state,  as  compared  to  the  correspond¬ 
ing  quantities  for  the  other  atoms. 

Comparison  of  the  computed  concentrations  in 
Fig.  6  and  those  measured  experimentally  [see 
Fig.  5(a)]  shows  that  a  kinetic  simulation  with  an 
assumed  anthracene-air  mole  ratio  of  5 1 1  to  10:1 
provides  a  reasonable  representation  of  the  experi¬ 
mental  results.  The  H  3  concentration  is  found  to 
be  the  largest,  as  was  observed  experimentally.  The 
computed  excited-state  concentrations  decay  at  a  rate 
similar  to  what  was  observed  experimentally.  The  ma¬ 
jor  difference  between  the  experimental  determined 
and  computed  concentrations  is  the  relative  concen¬ 
tration  of  the  C*5  state,  which  was  computed  to  be 
comparable  to  that  of  the  H*3  state,  while  the  C*5 
to  II  3  concentration  ratio  was  determined  experi¬ 
mentally  to  be  approximately  0.5  at  short  delays.  It 
is  possible  that  we  have  underestimated  the  uncer¬ 
tainty  in  experimentally  determined  the  C*5  relative 
concentration.  This  state  emits  in  the  UV,  far  from  the 
wavelengths  of  the  other  emission  lines,  which  occur 
in  the  red.  The  sensitivity  of  our  detection  system  is 
much  less  in  the  UV,  by  a  factor  of  ~40  at  248  nm 
(wavelength  of  the  C*5  line),  than  in  the  red,  at 
656  nm  (wavelength  of  the  H*3  line). 

C.  2,4-Dinitrotoluene 

Our  second  organic  sample  was  chosen  to  be  2,4- 
dinitrotoluene,  which  has  a  large  N  and  0  mole  frac¬ 
tion,  similar  to  that  in  many  explosives.  Figure  7(a) 
presents  the  determined  time  dependence  of  the  re¬ 
lative  concentrations  of  the  excited,  emitting  states 
of  the  light  atoms  and  A1  for  LIBS  of  residues  of 
the  compound  on  A1  foil.  A  laser  pulse  energy  of 
10  mJ  was  employed.  For  the  light  atoms,  the  deter¬ 
mined  concentrations  and  time  profiles  of  the  emit¬ 
ting  states  of  C,  N,  and  0  are  similar,  while  the  H*3 
concentration  is  approximately  a  factor  of  2  larger.  As 
we  observed  for  LIBS  of  bare  A1  foil  and  of  anthra¬ 
cene  residue,  the  A1  4S  concentration  decays  much 
more  slowly  and  persists  for  longer  delay  times. 

Figure  7(b)  presents  the  temporal  dependence  of 
the  plasma  temperature,  determined  from  the  ratio 
of  the  A1  394.40  and  308.22  nm  lines.  The  tempera¬ 
ture  was  found  to  be  somewhat  lower  than  that  found 
for  LIBS  of  anthracene.  It  should  be  noted  that  the 
laser  pulse  energy  was  10  mJ,  lower  than  that  em¬ 
ployed  for  LIBS  of  anthracene  (20  mJ). 

Kinetic  simulations  of  2,4-dinitrotoluene-air  mix¬ 
tures,  following  the  procedures  described  in  Section  2, 
were  carried  out.  As  was  done  for  simulations  of 
anthracene-air  mixtures,  A1  was  not  included  in 
the  reaction  mixture.  The  temperature  profile  in 


Fig.  7.  (Color  online)  (a)  Relative  populations  of  the  excited, 
emitting  levels  (in  the  notation  of  Table  1)  determined  from  the 
measured  intensities  (corrected  for  the  wavelength-dependent  de¬ 
tection  sensitivity)  of  atomic  emission  lines  for  the  355  nm  laser 
irradiation  (pulse  energy  10  mJ)  of  2,4-dinitrotoluene  residue  on 
A1  foil  in  ambient  air.  (b)  Plasma  temperature,  derived  from  the 
relative  intensities  of  the  A1  394.40  and  308.22  nm  lines,  as  a  func¬ 
tion  of  time  after  the  355  nm  laser  pulse  (10  mJ)  for  LIBS  of  2,4- 
dintirotoluene  residue  on  A1  foil. 

the  simulations  was  taken  from  the  fit  [smooth  line 
in  Fig.  7(b)]  to  the  temperatures  estimated  from  the 
A1  atomic  line  ratios.  As  in  Subsection  4.B,  we  took 
complete  atomization  and  partial  ionization  (see  the 
caption  of  Fig.  8  for  details)  as  our  assumption  for  the 
initial  concentrations.  With  this  assumption,  prompt 
formation  of  the  excited,  emitting  state  was  found  in 
the  simulations. 

Figure  8  presents  computed  time-dependent  con¬ 
centrations  of  the  excited,  emitting  light  atoms  for 
pure  2,4-dinitrotoluene  and  an  assumed  1:2.5  2,4- 
dinitrotoluene-air  mole  ratio.  The  H*3  state  is  com¬ 
puted  to  have  the  highest  concentration  in  both 
simulations,  while  the  second  most  populated  state 
is  predicted  to  be  the  C*5  state.  In  the  simulation 
of  pure  2,4-dinitroluene,  the  concentrations  of  the 
emitting  N  and  O  states  are  more  than  a  factor  of 
2  lower  than  that  of  H*3.  We  also  see  in  this  case  that 
the  O  atom  concentrations  are  predicted  to  be  ap¬ 
proximately  a  factor  of  2  larger  than  those  of  the 
N  atom  states.  This  is  consistent  with  the  larger  mole 
fraction  of  O  than  N  in  the  pure  compound.  The  si¬ 
mulations  of  the  2,4-dinitrotoluene-air  mixture  yield 
higher  relative  concentrations  of  the  excited  N  and  O 
atom  states,  as  expected.  In  this  case,  the  excited  N 
and  O  state  concentrations  are  comparable.  This  is  a 
result  of  the  greatly  different  N  to  O  mole  ratio  in  air 
than  in  2,4-dinitrotoluene.  Also  plotted  in  Fig.  8  is 
the  computed  electron  density.  As  in  Fig.  6,  the 
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Fig.  8.  (Color  online)  Plots  of  the  time-dependent  concentrations 
of  the  excited,  emitting  atomic  levels  (in  the  notation  of  Table  1)  for 
a  2,4-dinitrotoluene-air  mixture,  computed  with  the  kinetic  model 
described  in  Section  2  and  the  temperature  plotted  in  Fig.  7(b).  The 
assumed  initial  conditions  and  species  mole  ratios  were  as  follows: 
(a)  pure  2,4-dinitrotoluene,  with  assumed  complete  atomization 
and  10%  ionization  of  C  [C,  6.3;  C+,  0.7;  E,  0.7;  H,  6;  N,  2;  O, 
4],  and  (b)  2,4-dinitrotoluene-air  mixture  at  a  1:2.5  mole  ratio, 
with  assumed  complete  atomization  and  10%  ionization  of  N  [C, 
7;  H,  6;  N,  4.95;  N+,  0.55;  E,  0.55;  O,  6].  The  electron  concentration 
is  scaled  as  indicated  in  the  panels. 


electron  density  is  computed  initially  to  increase  and 
then  to  decay  more  slowly  than  the  excited-state  con¬ 
centrations. 

The  excited-state  concentrations  computed  with 
an  assumed  2,4-dinitrotoluene-air  mole  ratio  of 
1 : 2.5  [Fig.  8(a)]  provide  a  reasonable  representation 
of  the  experimentally  measured  excited-state  con¬ 
centrations  plotted  in  Fig.  7(a).  While  the  mole- 
cule-to-air  mole  ratio  in  the  optimum  simulations 
was  a  factor  of  ~2  different  for  anthracene  and 
2.4-dinitrotoluene,  the  atom  mole  ratios  are  similar. 
As  with  the  comparison  of  observed  and  computed 
concentrations  for  anthracene  [Figs.  5(a)  and  6(a),  re¬ 
spectively],  the  major  difference  is  the  relative  con¬ 
centration  of  the  C*5  state,  which  is  computed  to 
be  lower  than  is  measured  experimentally.  This  dis¬ 
crepancy  could  be  the  result  of  the  poor  UV  sensitiv¬ 
ity  of  our  detection  system,  and  the  consequent  error 
in  correcting  the  raw  intensities. 


5.  Conclusion 

In  this  paper,  we  have  presented  a  kinetic  model  of 
the  physical  and  chemical  processes  determining  the 


concentrations  of  the  excited,  emitting  atomic  states 
responsible  for  the  atomic  emission  signatures  in 
LIBS  spectra.  To  provide  a  complete  picture  of  the 
LIBS  process,  this  model  would  need  to  be  coupled 
with  a  theoretical  description  for  the  ablation  and 
the  hydrodynamics  of  the  plasma  plume.  It  is  also 
known  that  particles  are  often  ejected  from  surfaces 
in  laser  ablation,  and  this  process  would  need  to  be 
included  in  a  full  description  of  the  LIBS  process.  In 
this  work,  the  temperature  and  electron  density  as  a 
function  of  time  were  estimated  from  ratios  of  atomic 
line  intensities  and  Stark  broadening,  respectively.  It 
was  also  necessary  to  assume  the  initial  concentra¬ 
tions  of  the  various  species,  in  order  to  describe 
the  effect  of  the  laser  interaction  on  the  ablated 
material  and  the  ambient  atmosphere. 

We  have  compared  the  predictions  of  this  model 
with  recorded  LIBS  spectra  of  several  samples,  in¬ 
cluding  bare  A1  foil  and  two  organic  compounds,  an¬ 
thracene  and  2,4-dinitrotoluene,  in  ambient  air.  Our 
kinetic  simulations  predict  excited-state  concentra¬ 
tions  and  their  time  profiles  which  are  in  reasonable 
agreement  with  relative  concentrations  derived  from 
the  experimentally  measured  line  intensities  in  the 
LIBS  spectra.  We  also  were  able  to  obtain  an  esti¬ 
mate  of  the  degree  of  air  entrainment  into  the  plas¬ 
ma  plume. 

The  modeling  includes  several  components.  The 
first  is  the  matrix  of  rate  constants  and  their  depen¬ 
dence  upon  temperature,  and  in  some  case  their 
pressure  dependence.  Experimentally  determined  re¬ 
action  rate  constants  have  generally  been  determined 
for  temperatures  <3000  K,  and  these  rate  constants 
need  to  be  extrapolated  up  to  plasma  temperatures. 
Fortunately,  we  would  expect  chemical  compounds 
to  be  completely  atomized  at  these  high  temperatures. 
A  more  uncertain  component  of  the  model  is  the  need 
to  assume  initial  concentrations  to  describe  the  disso¬ 
ciation  and  ionization  caused  by  the  laser-matter 
interaction.  It  is  thus  gratifying  that  our  kinetic  mod¬ 
eling  does  yield  computed  relative  concentrations  in 
reasonable  agreement  with  the  experimentally  deter¬ 
mined  values.  Measurement  of  absolute  concentra¬ 
tions  would  provide  an  even  more  stringent  test  of 
the  kinetic  model;  unfortunately,  such  an  experiment 
is  fraught  with  significant  likely  systematic  errors. 
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